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Abstract 

In order to determine polarized parton distributions we have made a 
new NLO QCD fit using all experimental data on spin asymmetries 
measured in the deep inelastic scattering on different nucleon targets. 
The functional form of such densities is based on MRST2001 results for 
unpolarized ones. We get for polarization of quarks (at Q 2 = 1 GeV 2 ): 
An = 0.87, Ad = -0.38, As = -0.04. The total quark polarization 
is rather big and we obtain: AS = 0.44. As a result of our fit we 
get = gA = 1-24, the value which is close to experimental number. 
With negligible As and rather big AS (comparable to as) the results 
of our new fit are quite different in character from other fits. 
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The experimental data on deep inelastic scattering of polarized leptons 
on polarized nucleons was collected for many years in experiments made in 
SLAC [0], CERN § and DESY §. The results were analyzed by many 
groups and next to leading order (NLO) QCD polarized parton distributions 
were determined ||, [5|, || [7]]. In the present paper we do not consider any 
new experimental data. All the existing experimental data on polarized deep 
inelastic scattering were already included in our latest fits || fjj. However 
our method of determination of spin densities depends very strongly on the 
parton distributions for unpolarized case. We assume that the asymptotic 
behaviour of our polarized parton distributions is determined (up to the 
condition that the corresponding parton densities are integrable) by the fit 
to unpolarized data. 

Recently a new determination of unpolarized parton densities performed 
by Martin, Roberts, Stirling and Thorne (MRST2001) || was published. 
These distributions have substantially modified (compared to older MRST98 
U fit) small x behaviour of valence u quark and gluon, as well as gluon 
density is not positive for all values of x variable. We want to update our 
fit and to check how the new functional form of parton distributions influ- 
ences it. We will follow the method presented in |5|, 0] where we use MS 
renormalization scheme in QCD. Despite of the fact that we only modify 
the functional dependence of the fitted (at Q 2 = 1 GeV 2 ) parton densities 
the obtained polarizations of quarks (i.e., Am, Ad and As) and gluons (AG) 
are significantly changed and the value of total quark polarization AS is in- 
creased. We will use data for spin asymmetry at given x and different Q 2 
(431 experimental points). 

Experiments on unpolarized targets provide information on the quark 
densities q(x, Q 2 ) and G(x, Q 2 ) inside the nucleon. These densities can be 
expressed in term of q ± (x,Q 2 ) and G ± (x,Q 2 ), i.e. densities of quarks and 
gluons with helicity along or opposite to the helicity of the parent nucleon: 

q = q + + q-, G = G + + G~ . (1) 

The polarized parton densities, i.e. the differences of q + , q~ and G + , G~ 
are given by: 

Aq = q + -q-, AG = G + - G~. (2) 

We will try to determine Aq(x,Q 2 ) and AG(x,Q 2 ), keeping in mind re- 
lations between eqs. (1) and (2). 
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The formulas for unpolarized quark and gluon distributions determined 
(at Q 2 = 1 GeV 2 ) in the fit performed by Martin, Roberts, Stirling and 
Thorne [| (they use A^g 4 = 0.323 GeV and a s (M 2 ) = 0.119) are: 

u v (x) = 

d v (x) = 

2u(x) = 

2d(x) = 

2s(x) = 

G{x) = 

where: 

M(x) = 0.222x" L26 (l-x) 7 - 10 (l + 3.42v^+10.3x), 
6(x) = 1.195x a24 (l-x) 9 - 10 (l + 14.05x-45.52x 2 ). (4) 

We will split q and G, as was already discussed in [|, |5], [7| , into two parts 
in such a manner that the quark distributions q ± {x, Q 2 ) remain positive. Our 
polarised densities for valence quarks, sea antiquarks (the same distribution 
we take for sea quarks) and gluons (at Q 2 = 1 GeV 2 ) are parametrised as 
follows: 

Au v (x) 
Ad v (x) 
2Au(x) 
2Ad(x) 
2As(x) 
AG{x) 

where: 

AM(x) 
AM S 
A5(x) 
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0.158oT a75 (l - x) 3 - 33 (l + 5My/x + 55.49a;), 

0.040aT a73 (l - a;) 3 - 88 (l + 52.73^ + 30.65x), 

0.4M(x) - 6(x), (3) 

0AM(x) + 5(x), 

0.2M(x), 

1.90aT a91 (l - a;) 3 - 70 (l + l.26y/x - 1.43s) - 0.21x" L33 (l - a;) 10 . 



= x a75 (l - a;) 3 ' 33 (ai + a 2 y/x + a 4 x), 

= x - -™(l- x )™ 8 (b l + b 2 Vx~ + b 3 x), 

= 0AAM(x) - A5(x), (5) 

= 0.4AM(x) + A5(x), 

= 0.2AM s (x), 

= x-°- 91 (l - x) 3 - 70 (d 1 + d 2 ^/x~ + d 3 x) + x-° 83 (l - x) w d A . 



= x-°- 76 (l-x) 7 - 10 ( Cl + c 2V ^), 

= x- a76 (l-x) 7 - 10 (c ls + c 2sV ^), (6) 
= x a24 (l - x) 9 - 10 c 3 (l + 14.05x - 45.52x 2 ). 



In order to have finite polarization for partons we use less divergent distri- 
butions at x — > for sea quarks and gluons (in the second term). We also 
define: 

Am = Au v + 2Au, 

Ad = Ad v + 2Ad, (7) 
As = 2As, 

and 

AS = Au + Ad + As, 
a 8 = Au + Ad -2 As, (8) 
a 3 = Am — Ad = 

In order to determine the unknown parameters in the expressions for po- 
larized quark and gluon distributions (eqs.(5,6)) we calculate the spin asym- 
metries (starting from initial value Q 2 = 1 GeV 2 ) for measured values of Q 2 
and make a fit to the experimental data on spin asymmetries for proton, neu- 
tron and deuteron targets. The spin asymmetry A±(x, Q 2 ) can be expressed 
via the polarized structure function gi(x, Q 2 ) as: 

where R = [F 2 (l + 7 2 ) - 2xF 1 ]/2xF 1 , whereas 7 = 2Mx/Q (M stands for 
proton mass). We will take the new determined value of R from the [|1(] 



In calculating gi(x, Q 2 ) and F 2 (x,Q 2 ) in the next to leading order we use 
procedure described in |5|, []]. Having calculated the asymmetries according 
to eq.(9) for the value of Q 2 obtained in experiments we can make a fit to 
asymmetries on proton, neutron and deuteron targets. The value of CI3 is not 
constrained in such fit. We will also do not fix a 8 but we put a constraint 
on its value. Simply we will add it as an extra experimental point (from 
hyperon decays one has a§ = 0.58 ± 0.1, where we enhance (to 3a) an error). 

Not all parameters are important in our fits. We will assume that C\ s = C\ 
(i.e. the most singular terms for strange and non-strange sea contributions 
are equal). Such assumption practically does not change the value of x 2 but 
improves x 2 I^df- In this case we get the following values of parameters (at 
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Q 2 = 1 GeV 2 ) from the fit to all existing data for for spin asymmetries: 



at = 


0.16, 


a 2 = 


-2.45, 


a 4 = 


9.89, 


h = 


-0.04, 


b 2 = 


-1.03, 


h = 


0.65, 


Cl = 


-0.05, 


c 2 = 


2.28, 


c 3 = 


1.20, 


Cu = 


Cl, 


C2s = 


-0.33, 






di = 


8.67, 


d 2 = 


2.20, 


d 3 = 


-39.8. 



The quark and gluon distributions obtained from our fit lead for Q 2 = 
1 GeV 2 to the following integrated (over x) quantities: Am = 0.87, Ad = 
-0.38, As = -0.04, Au v = 0.77, Ad v = -0.67, 2Au = 0.10, 2Ad = 0.29. 

We have positively polarized sea for up and down quarks and small neg- 
atively polarized sea for strange quarks. One can see substantial breaking of 
SU(2) symmetry in a sea. We got the value of a 3 = g& = 1.24, the number 
which is very close to experimental figure (1.267 ± 0.003 [|TT|j). 

As is seen from eqs.(5) and (6) the small x behaviour of expressions for 
valence u, d quarks and sea contribution AM have very similar behaviour at 
small values of x. So one cannot expect that the splitting into valence and 
sea contributions will be well determined in the fit. The integrated values 
of quantities that are well determined i.e., Am = 0.87, Ad = —0.38, As = 
—0.04, AS = 0.44 and qa = 1-24 can be compared with the values from 
previous fit (called A' 2 in §) Am = 0.77, Ad = -0.57, As = -0.19, AS = 
0.01 and gA = 1.34. More singular behaviour of Au at small x in the present 
fit does not only directly influence the value of Am. We have some increase 
in Am and relatively big increase in Ad and As resulting in the higher value 
of AS. It seems that improvement of the MRST2001 || in comparison to 
MRST98 changes significantly character of obtained results. On the other 
hand the value of \ 2 — 376.7 for the present fit is higher then in our previous 
one 0, despite of the fact that we have used the same experimental data. 

The results in the "measured" region (i.e., for 0.003 < x < 1) are: Am = 
0.76, Ad = -0.32, As = -0.03, AS = 0.41 and g A = 1-08 which should be 
compared with our previous results of fit A' 2 Au = 0.80, Ad = —0.43, As = 
—0.11, AS = 0.26 and gA = 1.23. In present fit we have the value of AG = 
25.7 integrated in the whole region of x and AG = —1.8 in the "measured" 
region of x, which have to be compared with results of a A' 2 fit: AG = —0.19 
in the whole region and the same value in the region: 0.003 < x < 1. 

As we already mentioned before the assumption that AG should be inte- 
grable was used by us in basic fit and the most singular part of G was removed 
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from AG parametrization. We have not assumed positivity for gluon con- 
tribution because already in unpolarized case G{x) was negative for small x 
values. When we take for AG(x) the same functional form as for G(x) leav- 
ing most singular term for small x (not integrable) we get fit with smaller % 2 , 
namely 372.8 (there are some changes in parameters but the character of the 
fit does not change). We get: Au = 0.87, Ad = -0.38, As = —0.04, AS = 
0.45 and gA = 1.24. The value of gluon polarization (AG) is —0.79 in the 
measured region and of course infinite in the whole region. We have made also 
a fit where we have not assumed any positivity conditions for parton densi- 
ties. It means that functional form for such densities was used with arbitrary 
(unconstrained) parameters in eqs.(5) and (6). We get, as expected, much 
smaller \ 2 value, namely 360.1 (smaller than in our basic fit). The integration 
in the whole x region gives: Am = 0.84, Ad = -0.38, As = -0.04, AS = 0.41 
and gA = 1.22. AG = —12.2. When we neglect the second term in 
gluon distribution (i.e., when we put d 4 = 0) one gets \ 2 — 384.2 and: 
An = 0.85, Ad = -0.40, As = -0.06, AS = 0.39, AG = -14. If we 
put (for Q 2 = 1 GeV 2 ) gluon polarization arbitrarily equal to zero (i.e. 
d\ = d 2 = d% = d± = 0) we get substantial increase in x 2 = 393.2, whereas: 
Am = 0.84, Ad = -0.42, As = -0.07, AS = 0.36 and g A = 1.26. 

It seems that the values of Au, Ad, As, AS and qa do not change much 
when we use different assumptions about gluon contribution. We get in our 
fits small value of As and AS comparable with as (such conclusions remind 
Ellis- Jaffe results |12|). This is what differs our present results (coming from 
updated MRST2001 {§ fit) from our previous fits §, §. 

Our fitting procedure also determines + and — helicity components of 
parton densities. These components and the total polarized distributions 
for quarks of different flavours and gluons (calculated at Q 2 = 1 GeV 2 ) are 
presented in Fig.l. 

In Fig. 2 the values of functions of Au(x) , Ad(x) , AT,(x) and AG(x) ob- 
tained in the present fit are compared with the corresponding ones from our 
previous fit called A' 2 0. One can see significant changes (except maybe Au) 
in these distributions. 

In Fig. 3 we present the comparison of structure functions g\, gf and g™ 
(calculated at the values of Q 2 corresponding to data) with the correspond- 
ing experimental points. For comparison, the curves for polarized structure 
functions obtained in || from fit A' 2 are also presented. The increases of g\ 
(however within experimental errors) are seen only for very small x values 
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Figure 1: The quark (for different flavours) and gluon densities versus x ob- 
tained from our fit. The solid lines represent parton polarization, whereas 
dashed (dotted) lines correspond to + (—) helicity component of parton dis- 
tribution. 
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Figure 2: The parton densities versus x obtained from our fit (solid lines) 
compared with the results from our previous fit J^/ (dashed lines). 
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Figure 3: The comparison of our predictions for gi(x,Q 2 ) versus x with 
the experimental data from different experiments. Solid curves are calculated 
using the distributions from our fit, the dashed ones are calculated using the 
parameters of fit A 2 from j^. 
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for data from SMC experiments (on proton and deutrium targets) and for 
El 54 data on neutron target. 

We have made fits to precise data on spin asymmetries on proton, neutron 
and deuteron targets. Our model for polarized parton distributions is based 
on new MRST2001 fit || to unpolarized data. Because of different functional 
form of fitted parton densites (e.g. more singular behaviour at x = for u 
quark) the new fit is different in character from the other fits 0, [5], [|. We have 
got relatively small value of As and relatively high value of AE (comparable 
with a§) with gA (unconstrained) very close to experimental value. Gluon 
contribution comes out relatively high and is (as usual) not very reliable. It 
seems that in our model the problem of so called "spin crisis" practically 
disappears. 
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